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Abstract 

We apply the method of slowly- varying amplitudes of the electrical 
and magnet fields to integro-differential system of nonlinear Maxwell 
equations. The equations are reduced to system of differential Non- 
linear Maxwell amplitude Equations (NME). Different orders of dis- 
persion of the linear and nonlinear susceptibility can be estimated. 
This method allow us to investigate also the optical pulses with time 
duration equal or shorter to the relaxation time of the media. The elec- 
tric and magnetic fields are presented as sums of circular and linear 
components. Thus, NME is written as a set of Nonlinear Dirac Equa- 
tions (NDE). Exact solutions of NDE with classical orbital momentum 
£ = 1 and opposite directions of the spin (opposite charge) j = ±1/2 
are obtained. The possible generalization of NME to higher number 
of optical components and higher number of £ and j is discussed. Two 
kind of Kerr type media: with and without linear dispersion of the 
electric and the magnet susceptibility are consider. The vortex solu- 
tions in case of media with dispersion admit finite energy while the 
solutions in case of a media without dispersion admit infinite energy. 
PACS 42.81.Dp;05.45Yv; 42.65.Tg 

1 Introduction 

In the last four decades there has been considerable interest in the 
nonlinear generalizations of the quantum field equations IS and 
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in the possibility of obtaining exact stationary solitary solutions of 
the field equations As a rule, the nonlinearity has been intro- 

duced in an ad-hoc fashion in the Klein-Gordon equation and also for 
all four spinor components of the Dirac equations. For the usual case 
of a cubic nonlinearity, exact localized solutions are not found. Our 
present work, reported in this paper, shows that the optical analogy 
of nonlinear Dirac equations leads to a nonlinear part only the first 
coupled equation of Nonlinear Dirac Equations. This result allows to 
solve the NDE by separation of variables and to obtain solutions repre- 
senting optical vortices with classical momenta one and spin one-half. 
The initial investigation of optical vortices began with a scalar theory, 
based on the well-known 2D+1 paraxial Nonlinear Schrodinger equa- 
tion (NSE) IHl [Tj . The existence of optical vortices was predicted in 
the self-focusing regime, but as it was shown in many papers, the so- 
lutions obtained are modulationally unstable. In spite of this, various 
interactions (attraction, repulsion, fusion) have been observed. The 
scalar paraxial approximation is valid for slowly varying amplitudes 
of the electrical field in weakly nonlinear media. As it was pointed out 
in Refs.|Hlin| this theory is not valid for a very intense narrow pulses. 
The first generalization of the scalar paraxial theory of optical vortices 
is based on investigation of so called spatio-temporal scalar evolution 
equations ^1 ^2 . For all cited theories no exact solutions have 
been found, but numerical and energy momentum techniques are used. 
The existence of exact stable vortex solutions of these types of non- 
linear equations was finally discovered with the vector generalization 
of the 3D+1 Nonlinear Schrodinger equation jj^. It also have been 
shown numerically that these vortices are stable at distances compa- 
rable to those where localized solutions of the one component scalar 
equation self-focusing rapidly. All of these exact vortices are a com- 
bination of linearly-polarized components and have spin i = 1. To 
extend the theory to vortices with spin j = 1/2 we return to an anal- 
ogy between the Maxwell and Dirac equations . As it is shown in Ref. 
jl4j this analogy is possible only if the electrical and magnet com- 
ponents are represented as a sum of linear and circularly polarized 
components. 
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2 Maxwell's equations with non-stationary 
linear and nonlinear polarization 

Consider the Maxwell's equations in the next two cases: 

1. A source-free medium with non-stationary linear and nonlin- 
ear electric polarization and non-stationary magnetic polarization (the 
case with dispersion). 

2. A source- free medium with stationary linear and nonlinear elec- 
tric polarization and stationary magnetic polarization (the case with- 
out dispersion). 

For these cases, the Maxwell's equations can be written: 



VxE = — , (1) 



c 



dt ' 



V-D = 0, (3) 

V ■ B = V ■ H = 0, (4) 

D = P'^ + 4TrPne. (5) 
The linear magnetic polarization corresponds is [18,19]: 



B = H + 4TrMun, (6) 

where E and H are the electric and magnetic intensity fields, D and B 
are the electric and magnetic induction fields, p'*"^, P^^ are the linear 
and nonlinear polarization of the medium respectively and is the 
linear magnetic polarization. The magnetic polarization (magnetiza- 
tion) Miin is written as the product of the linear magnetic suscepti- 
bility r]^^^ and the magnetic field H. The nonstationary linear electric 
polarization can be written as: 
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p'-= J (^5{t-T)+A7rx^'Ht-T))E{T,x,y,z)dT 

— oo 

t 

= j eoit-T)E{T,x,y,z)dT, (7) 

— oo 

where x^^^ and Eq are the hnear electric susceptibihty and the dielectric 
constant respectively. Similar expression describes the dependence of 
B on H in the case of nonstationary linear magnetic polarization j^S] : 



t 



B = J (5{t-T)+ATTr]^^^ {t-T))H{T,x,y,z)dT 



-oo 

t 



Ho{t-T)H{T,x,y,z)dT, (8) 



where rj^^^ and /io are the linear magnetic susceptibility and magnetic 
permeability respectively. The magnetic susceptibility of the main 
part of the dielectrics ranges from 10~^ — 10~^ and usually decreases 
with the increasing of the frequency. In the following, we will study 
such media with nonstationary cubic nonlinear polarization, where 
the nonlinear polarization in the case of one carrying frequency can 
be expressed as: 

— oo — oo — oo 

xE{ti,x, y, z)E*{t2,x, y, z)E{t3,x, y, z)dTidT2dT^. (9) 
The causality request the next conditions on the response functions: 

e(t-r) = 0; x^^^ - n, t - ra, i - rg) = 0, 

t - T < 0; t - Ti < 0; i = 1, 2, 3. (10) 

That why we can prolonged the upper integral boundary to infinity 
and to use standard Fourier presentation |18j : 
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J Eq^t — t) exp {iuJT)dT = J £o{t — t) exp {iuJT)dT (11) 

-oo — oo 

t t t 

III X^^^ (*-^i'*-^2,t-r3)(iri(ir2dT3 = 

— oo — oo — oo 

+00 +00 +00 

III ^^^^ (*-'^i'*-^2,t-T3)(iridr2dT3. (12) 



-00 —00 —00 



The spectral presentation of linear optical susceptibility £q{uj) is 
connected to the nonstationary optical response function by the next 
Fourier transform: 



+00 



£q{uj) exp {iojt) = j £o{t — t) exp (iujT)dT, (13) 
—00 

Similar expression for the spectral presentation of the non-stationary 
magnetic response ^0(1^): 



fiQ{uj) exp {—iujt) = J fj,Q{t — t) exp {—iuJT)dT . (14) 
—00 

and nonlinear optical susceptibility x^^^ 



X 



+00 +00 +00 

(3) 



III x'^^Ht-ri,t-T2,t-T3) 



-00 —00 —00 



X exp (i (cj (ti + T2 + T^)))dTidT2dT-i (15) 

can be written. It is important to point here the remark of Akhmanov 
at all in |^ , that such nonstationary representation applied to slowly 
varying amplitudes of electrical and magnet fields is valid as well when 
the optical pulse duration of the pulses Iq is greater than the charac- 
teristic response time of the media tq (to >> ^o), as when the time 
duration of the pulses are equal or less than the time response of the 
media (to < ''"o)- We will discuss this possibility in the process of 
deriving of the amplitude equations. 
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3 Deriving of the amplitude equations 



In this section we derive the slowly varying amplitude approximation 
in the standard way, as it was done in J7j, JH]. For the case of 
Maxwell equations with linear and nonlinear dispersion Q-® we 
define the electric and magnetic field amplitudes with the relations: 



fields, the optical frequency and the real spatial phase functions re- 
spectively. In case of monochromatic and quasi-monochromatic fields 
the Stokes parameters can be constructed from transverse components 
of the wave field [W\ I2()j . This leads to two component vector fields in 
a plane, transverse to the direction of propagation. For electromag- 
netic fields with spectral bandwidth (our case) the two dimensional 
coherency tensor cannot be used and the Stokes parameters cannot 
be found directly. As it was shown by T.Carozzi and all in |2nj . 
using high order of symmetry (SU(3)), in this case six independent 
Stokes parameters can be found. This corresponds to three compo- 
nent vector field. Here is investigated this case. The increasing of 
the spectral bandwidth of the vector wave, increases also the depo- 
larization term (component, normal to the standard Stokes coherent 
polarization plane). 

We use the Fourier representation of the response functions (| 1, 
H15|l and of the amplitude functions A (and C) to obtain next ex- 
pressions for the first derivatives in time of the linear polarization, 
nonlinear polarization and magnetic induction fields: 




(16) 





iexp (i {uJot - g{x,y,z))) 



c 



dt 




(18) 



— oo 
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X / \A[x,y,z,LU - LUo)\ A[x,y,z,LU - loo) 

—oo 

X exp {i{Lo — LOo)t) du. (19) 

1 dB {x,y,z,t) 



X 
—oo 



c dt 

+ 00 

ufio (w) 



iexp{i {uJot-q{x,y, z))) 



/(jj ioj I 
C{x, y,z,LO — uJo) exp — uJo)t) du. (20) 



At this point we restrict the spectrum of the ampHtude of electrical 
and magnet fields by writing the wave vectors fci,ni,2 in a Taylor series: 



ujip (u;) 
ki [Lo) = = ki {cjo) + 



ki (wo) + — (a; - a;o) + l-kl {u - uof + (21) 



Sttc^X^'^ (^) 

^nZ (t^) = — 



,0/ ^. 9{knl (wo)) . s . 



k°nii'^o) + — i'^-u;o) + ...., (22) 

Vnl 



k2 [Lo) = — = k'i (wo) 
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^2 ('^o) + — (o; - (Jo) + ^^2 - '^o)^ + (23) 

where Vi and k^;i = l,n/,2 have dimensions of group velocity, and 
nonlinear addition to the group velocity and dispersion respectively. 
The nonlinear wave vector is expressed: 



Ki = SvrwoX^ n^^o)/c = 2 = kin2, (24) 

c 

where 



n2(uJo) = , (25) 

^0 

is the nonlinear refractive index. It is important to note here that 
we do not use any approximation of the response function. There is 
only one requirement of the spectral presentations (fT5|) . (|13|) and (|Ti|l 
of the response functions: to admit first and second order derivatives 
in respect to frequency ( to be smooth functions). The restriction is 
only in respect of the relation between the main frequency ujq and time 
duration of the envelope functions to determinate from the relations 
(PT|) , and (conditions for slowly varying amplitudes) . Putting 
eqn. (ETJ in (UHj), ^ in ^ and ^ in and keeping in mind 
the expressions for time derivatives of the spectral presentation of 
the amplitude functions, the first derivatives ((TH|l . (fT^ and (^0]) are 
presented in form: 

1 gP'^" {x,y,z,t) _(.,ox^ ^ ikl d^A\ 

c Ft -y^'^^n^ 

xexp{i{u;ot - g{x,y,z))), (26) 
47r dP^^'^ix,y,z,t) _ ( , 1 d\A\^A \ 

al ~ lKin2\A\ A-\ 

C dt \ Vnl Ot J 

xexp{i{uJot - g{x,y,z))), (27) 
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ldB{x,y,z,t) _ ( ro^ 1 dC ikld^C\ 

~ c Ft " y''^^ + ~ ) 

xeyip{i{ujot- q{x,y,z))). (28) 

Finally, from the Maxwell equations (0)-®, using eq. (|26|). eq. (|27|) 
and ea. H28() . and using in fact that V • ~ V • ~ |13j . we obtain 
the next system of Nonlinear Maxwell vector amplitude Equations 
(NME): 



VxA = <C--^-^^. (29) 



V X C = ik\A + —^-^^ + ik\n2 [A- A*] A 
vi at 2 at^ V / 



+ [ — + ^ ^ , (30) 



.til / dt 

V-A = 0, (31) 

V • C = 0, (32) 

if the gradient of the spatial phase functions g and q satisfied the 
relations: 

V5xi'=0, (33) 

and 

Vg X C = 0. (34) 

The phase functions whichever satisfied (|33|) - (|34|) are determinate in 
Section (10). 
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We investigate the case, when our vector fields are presented as a 
sum of circular and linear polarizing components: 



A — Aiin + Acir, (35) 
C = Clin + Ccir- (36) 

The nonlinear polarization admit different nonlinear refractive indexes 
in the case of linear and circular polarization \'21\ ( 71,2" 7^ '^2^). We 
will include this difference in the our rescaled equations, defining the 
next rescaled dependant variables: 

A = Aq^'Ah^ + ^ciri (37) 
C = C'o(C/m + C'cir); (38) 

{A't? = %iA^n', (39) 

77-2 

and independent variables; 

X = rox';y = roy' ; z = roz ,t = tot . (40) 
In additions the next constants are determinate: 

ai = k^ro; Pi = /cj ro/2tg; 71 = rohnf'' \Ao\'^ ; 

72 = nf l^ol'; 73 = t'inrl^olVc; 
6 = vi/v2 < 1; wi ~ ro/to; z = 1,2. 

The NME H29 |) -()32 |) in rescaled variables are transformed into the fol- 
lowing (the primes have been omitted for clarity): 



V xA = ia2C-6—-i(32^, (41) 



V X C 



. dA 
ot 



d'^A 



+ (72 + 73, 



dt 



(42) 
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v-i*=o, 



(43) 



V • (7 = 0. (44) 

We consider the case of slowly varying amplitudes approximation when 
the nonlinear constant is 71 = 1. Then the constants in fronts of last 
therm on right hand are 72 ~ 73 ~ 10~^ — 10^^. While the constants 
ai have typical values a 10^ — 10^ (aj ~ rofe,), the dispersion con- 
stants Pi have very small values j3i ~ 10~^ — 10~^ for picosecond and 
sub-picosecond pulses in the transparency region of nonlinear optical 
media. Neglecting the small dispersion terms [3i « 1 and the last 
nonlinear term for witch 72 ~ 73 << 1, the NME system H41() - ()44|) can 
be rewritten: 



V X A = ia2C - 6—, (45) 

yxC = iaiA + — +i-ii{A-A*]A, (46) 
at ^ 

V • i* = 0, (47) 

V • C = 0. (48) 



4 Dirac representation of NME 

To solve the NME we apply the method of separation of 

variables. The slowly varying amplitude vector of the electric field A 
and the magnetic field C are represented as: 

A{x,y,z,t) = F{x, y, z) exp (iAat) , (49) 
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C (x, y, z,t) = G (x, y, z) exp (iAat). (50) 
Substituting these forms into the NME (|45jl -(|48 jl we obtain: 

V X F = -iu2G, (51) 

V X G = iviF + i7i {f ■ F*) F, (52) 

V • F = 0, (53) 

V • G = 0, (54) 



where 1^1 = 01 + Aa; 1^2 = 5Aa — 02 > 0. In a Cartesian coordinate 
system, the vector equations are reduced to a scalar system 

of eight nonhnear wave equations. When the electric and magnetic 
fields are represented as a sum of a linear polarization component and 
a circular polarized one it is possible to reduce eqns. (|?T|) - H54() to a 
system of four nonlinear equations. Substituting Jl]: 



V'l = iFi = iF„ (55) 

i^2 = Fc = iF^ - Fy, (56) 

^3 = Q = -G„ (57) 

= Ge = -G, - iGj,, (58) 



into the nonlinear system H51|) - (|54|) . we obtain a stationary nonlinear 
Dirac system of equations (NDE): 
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d . a \ d ^ 

— + I— ^3 - — ^'4 = -I 

ox ay J oz 



(^i+^iEl^il') *2, (60) 



- ^TT ) ^2 + TT^i = -^^^2^-3, (61) 
ax ay / Oz 



d d\ d 

+ ^1 - ^2 = -W2^A. (62) 
ax ay / Oz 

This substitution allows to reduce the system of eight equations ((SJ)- 
(|54|) to a system of four scalar complex equations (|59jl - (|62jl . The 
system H59() - ()62|) is the optical analog of the nonlinear Dirac equations 
(NDE) . Note that the optical NDE are significantly different from the 
NDE in the field theory. The nonlinear part appears only in the first 
two coupled equations of the system. 



5 Vortex solutions with orbital mo- 
mentum 1=1 and spin j=l/2 

The symmetries of the NDE (|59|) - (I62|) are used to obtain exact vortex 
solutions. The NDE 1)59(1 - 1)621) have both spherical and spinor symme- 
try only in the case where the nonlinear part does not manifest the 

2 2 

angular dependence on the radial variable = F{r). 

i=l 

This type of solution can be found using with the following technique. 
Using Pauli matrices, we write the NDE system ()59() - ()62|) as: 



^1^1 + 71^ l^il^^ (63) 



a-P)T] = -iu2<j), (64) 
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where 

■"'O l\ f -i\ ( I 



1 / ' \ i / ' \ -1 



are the PauU matrices, P = (^^i ^) 



is the differential operator and 

are the corresponding spinors. After substituting eqn. (|64j) into 
ean.(|S5|l we obtain: 



-V2 {^1 +liYl^ fJ- (65) 



When there is no external electric or magnetic field, the operator on 
the left-hand side of eqn. (|65() is is the Laplacian operator A: 

p2 = . (a . p) = A. (66) 
From and we obtain: 

2 

1^21^11] + iy2ji^\T]i\'^r] + At] = 0. (67) 

1=1 

In case of spherical representation of the spinor equations (|65|) . there 
are two possibilities, I = and / = 1, that will permit only a radial 
dependence of the nonlinear part: 

2 

^|7?,(r,e,(^)|2 = P(r). (68) 
1=1 

For the I = case, the equations (|H7)l are transformed to equations 
with radial parts: 



^ 2 d^V "^^rj 
'^2i^iV + '^2li}^\Vi\ V+Q^ + --Q^ = 0. (69) 

i=l 

The scalar variant of these equations has been investigated in many 
papers but exact localized solutions have not been found. In the case 
/ = 1 we look for spinors in the next two forms: 



/ 7?(r)cos(0) \ 
\ f] (r) sin {6) exp {iip) j 
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and 



fj (r) sin (6) exp {—iip) 
—fj (r) cos {9) 



(71) 



As it will be seen later this corresponds to two opposite directions of 
the own orbital momentum j = ±1/2 (opposite charge). After substi- 
tuting solutions (|7()j) - (|71j) into equations (jlTTj) the following equation 
describing the radial dependence is obtained: 



,~,2~ d'^V 2dn 2 

^2^17] + U2li\v\ 1 + + = 0. 

or^ r or 



(72) 



The angular parts are the standard spherical harmonics with / = 1. 
This system has exact vortex de Broglie soliton solutions in the 
form: 



fj (r) 



a/2 exp {i^/uTu2'> 



(73) 



I r 

if z^27i = 1- The complete solutions for these two cases are written: 



V ^ -P(-xf^-) sin (g) exp (^y) 



(74) 



and 



IU2 \ 



'u\V2 exp {i^i'\V2r) exp i(^/i7[u^r) 



+ 



V 



-V2 exp(i^uiU2r) 



COS (0) 



(75) 



(76) 



\ iu2 \ r J 



(77) 
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There is another, more direct way for separating the variables of the 
spinor equations (|63() - H64|) . To illustrate this, we represent the NDE 
(|59 |) -(|62 |) in spherical variables: 



, . X / . cos 6 d i d \ ^ 

exp {-tip) sm9— + — ^4 

V or r oa r sm 0(p / 

+ ( ^1 - ^3 = f .1 + 71 E 1^.1^1 ^1, (78) 



i=l 



cos ( 



/ d cos 9 d i d \ 

exp {i(f) sm9— H — H ^tt- ^3 

V or r 09 r sm 9 Oip / 

a sin 61 a \ / A, ,2\ 

^4 = -^Ui+7iEl*irP2, (79) 



dr r 89 



i=l 



, . . f . n d cos9 d i d \ ^ 

exp {-1(f) sm9— H — ^tt- ^2 

V or r 09 r sm 9 0(p / 

/ sin9 d \ ^ ^ 

+ fcos^-----j^i = -zz.2^3, (80) 



, , , „ d cos 9 d i d \ ^ 

exp {lip) sm9— + — + —^jr- ^'i 

' or r o9 r sm 9 dip J 

d sin 6* a \ . ^ 

cos^- — U-a = (81) 

Or r 09 J 

We make the following two ansatzes for solutions to the system of 
nonlinear equations (|78|) - (|81|) : 



xl/^ = a (r) cos {9) 
^2 = a{r)sm{9)e'^ 

^2, = -ib{r) ^ ' 



= 0. 

and 



■^i = a (r) sin {9) e'^P 
^2 = -a (r) cos {9) 

*3 = ^^^^ 

\I/4 = -ib{r) . 
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Put or l|S5|) in equations ((75|) - (|ST|) we separate the variables. The 
following system of equations describing the radial dependence of the 
amplitudes are obtained: 

9a (r) 2 , , , , , , 

-^ + -a{r) = -U2h{r), 84 
or r 

-^=j,^a{r)+-ii\a{r)\'a{r). (85) 

In the system of equations ()84() - H85|) the nonlinearity appears only in 
the radial part, while for the angular part we have standard spherical 
spinors with spin j = ±1. Solving eans.H84|) and 1)851) we straightfor- 
wardly show that when V2I = 1 these equations admit the localized 
de Broglie solitons of eqns. H74() - ((75|l . 



6 Hamiltonian representation of the 
NDE. First integrals for vortex solu- 
tions with spin j = ±1/2 

It is not difficult to show that for the NDE system of eqns. H59|l- H62|l 
a Hamiltonian of the form: 

2 

H=(&-P)+Y.\^i\\ (86) 

i=l 

can be written. Using this, eqns. (fS^ -(jH ^ can be rewritten in the 
form: 

= e-^, (87) 

where e = {—ivi,—iui,—iu2,—ii'2) is the energy operator. Here we 
investigate the case where the nonlinear part of the equation is rep- 
resented as a number of spinors with a scalar sum that depends only 
on the radial component. 

2 

Y,\^,{T,e,^)\^ = F{r). (88) 

i=l 

We also introduce here the well known orbital momentum operator L, 
own orbital (spin) momentum 5, and the full momentum J, as well 
as: 
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L = fxP 



(89) 




(90) 



(91) 

It straightforward to show that the Hamiltonian l|86|) of eqn. (|87|) 
commutes with the operators and Jz (the z- projections must be x 
or y). Using these symmetries and the condition that the nonhnearity 
is of Kerr type, we can solve the NDE equations (|S7|) by a separation 
of variables technique. We look for solutions in the form: 



^i = a (r) QjiM (92) 
^2 = a (r) QjiM (93) 
M'3 = i6(r)17^,,^, (94) 

^4 = i6(0^^,7'M' (95) 

where i^jim is the spherical spinor, I + l' = 1, and a (r) and b (r) are 
arbitrary radial functions. Using the symmetries of ()87|1 and the fact, 
that the nonlinear parts depend on r, we separate variables and obtain 
the following system of equations for the radial part: 



5a(r) 1+X . ^ , . ^ .„„^ 

\ a{r) = -U2b{r) (96) 

or r 



^ + L^b (r) = uia (r) + ^\a {r)\''a (r) , (97) 
or r 

where 

X = /(/ + l)-i(j + l)-l/4. (98) 

Exclude 6(r) from the system (|^ - H97|) . we obtained the next equation 
for a(r): 
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d'^a 2da (1 + x)x , , ,2 n /nn\ 
^i^2a{r) + 0^ + -^ -2 « + ^2l\a\ a = 0. (99) 

Formally this equation admit exact "de Broglie" soliton solutions 
for arbitrary number of x- But as we remember our solutions are 
limited by the conditions (|88() the nonlinear part to depends only on 
the radial components. The condition (|88)) for a number x ^ 1 can be 
fulfilled also for a higher number of field on different frequencies. This 
case including also the parametric processes. The case of one carrying 
frequencies correspond to localized solutions with x = 1 ^-iid angular 
components I = 1 and j = ±1/2. In this case the system (|^H|) - (|^7|) 
becomes: 



'^!^SLl + _a{r) = -V2h{r) (100) 
or r 



^^ = ^^a{r)+^\a{r)\\{r). (101) 
or 

As was shown above, this system has exact radial solutions of the 
form (f7n|)- (|7^ : 



7 Experimental conditions 

There are some differences between the nonlinear conditions for local- 
ized solutions of the Vector Nonlinear Schrodinger Equation (VNLS) 
and localized conditions for the Nonlinear Maxwell Equations (NME). 
The nonlinear parameter for the VNLS is written |13) : 

Ivnls = ^0^0^2|^o|^ = 1- (102) 

For localized solutions in optical region the constant ranges from: 

= klrl K, 10^ - 10^, (103) 

which corresponds to a required nonlinear refractive index change of 
order of: 

nzlvlol^ ~ 10"^ - 10"^ (104) 
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On the other hand, the nonhnear condition for locaHzed solutions of 
the normahzed NDE (gnil-llHl) is: 

INME = ~ (5Aa - a2)kirQn2\AQ\^ = 1. (105) 

For the typical value of the constant 1^2 ~ 1 the next required nonlinear 
refractive index change appears: 

nsl^ol^ ~ 10"^ - 10"^ (106) 

Another difference between these two cases is that the solutions of the 
VNLS are comprised of linearly polarized components and that the 
dispersion of the nonlinear medium plays an important role. This leads 
to the fact that vortex solutions of the VNLS may be observed only in 
special dispersion regions of nonlinear media. The vortices of the NDE 
are a combination of linear and circular polarization components, and 
do not have this marked dependence on the dispersion and may be 
observed in the transparency region of a nonlinear media. 

8 Vortex solutions in a nonlinear Kerr 
type media without dispersion 

The case without linear dispersion of the electrical and magnetic sus- 
ceptibility will correspond to x^^^ = const and ry^^^ = const . We 
suppose now that the amplitude functions do not depend from time 
and will look for 3D+1 monochromatic electric and magnet field of 
kind: 

E{x, y, z, t) = M{x, y, z) exp {iioot), (107) 



i?(x, y, t) = N{x^ y, z) exp (iujot), (108) 

where M N and ujq are the amplitudes of the electric, magnetic fields 
and the optical frequency respectively. 

Substituting the relations H1U7() - H1U7() one obtain the next ampli- 
tude equations in the case of dispersionless media: 

V X M = ia2N, (109) 
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V X N = iaiM + j {^M ■ M*^ M, (110) 
V-M = (111) 

V-iV = (112) 

Again the electric and magnetic fields are represented as a sum of 
a linear polarization component and a circular polarized one. Substi- 
tuting: 

^1 = iMi = iM^ (113) 
11,2 = Mc = iM^ - My (114) 

^Ij^ = Ni = -iN, (115) 

'iP^ = N^ = -N^-iNy, (116) 

into the nonlinear system H109|) - ()112() . we obtain the same kind of 
stationary nonlinear Dirac system of equations (NDE) as (|59|) -(|62 |) : 



ox oy J oz 



(118) 



dx~'^dy) dz (119) 

but with different coefficients: = = C = ^^"f ■ 

Naturally, the NDE (fTT7 |) -(fnn |) admit the same kind of solutions (ffUl - 

m 
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9 Conditions for finiteness of the en- 
ergy of the vortex solutions 



In this section we will investigate the localized energy in the both 
cases; with and without of linear dispersion. This corresponds to two 
kind of optical vortices; with and without spectral bandwidth. In a 
dielectric media without dispersion the expression for the linear part 
of energy density is: 



E 



H 



[121] 



where e and ^ are constants. Substituting the vortex solutions H74|) - 
(|75|) in (|12H) and integrating over 3D space one obtain that quasi- 
monochromatic vortices in a dielectric media without dispersion admit 
infinite energy. Now we come to the case of slowly varying amplitude 
approximation and to calculation of energy of the vortex solutions 
of a media with linear electric and magnet dispersion. To prove the 
finiteness of energy of the vortex solutions ()74() - H75|) we start with the 
equations for averaged in time balance of energy density of electrical 
and magnet field [T3] : 



/dW\ 1 ( ^ dD* dD ^ dB* dB\ , , 

where D = Pun + AT^PnUn is a sum of the linear induction and the 
nonlinear induction of the electrical field. The calculations of the av- 
eraged energy of the optical waves in a dispersive media are worked 
out considering the first order of slowly varying amplitude approxi- 
mation of electrical induction (the same as in the NDE). The result 
comes to the old result of Brillouin 1921 for energy density of electrical 
field: 



1 (diujio) ^2 ^ d{ujfiQ) 



du> 



+ 



C 



(123) 



The conditions of electric constant: 



io > 0; > 

OLO 
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(124) 



is fulfilled in most of the dielectrics. The conditions of the magnetic 
constant can be present in two cases. In the first one: 



Ao>0;^>0. (125) 

This condition corresponds to the case, when there are not included 
any magnet relaxation processes. To include these processes is possi- 
ble only for ultrashort optical pulses with time duration of dimension 
of local structure of one paramagnetic media. The real part of mag- 
net susceptibility in this case is symmetric function in respect to the 
paramagnetic resonance and the derivative in respect to the frequency 
is always negative. 



Ao > 0; ^ < 0. (126) 
Under appropriate conditions this leads to: 



Ao>0;^i^^<0. (127) 

We can find spectral region and dispersion parameters when the 
next condition will be satisfied: 



d (c^Ao) d jujio) 

Thus, the linear (infinity) part of energy density is zero. The 
nonlinear part of averaged energy density is expressed in [21] and for 
the vortex solutions ()74() -H75 |) becomes: 



(WnHn) = n^l^f^i^^) + u;o7^|^|'3f?(^'). (129) 

otoo 

These results give the conditions for the finiteness of energy of the 
vortices. Integrating Wniin in the 3D space we obtain a finite value 
proportional to the main frequency u)q. 
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10 Spatial phase functions, Poynting 
vector and flow of energy 

The kind of the phase functions which satisfied (jHSJ-llSll is determi- 
nate for vortex solutions ()82jl with spin j = 1/2. Using again the 
relations between the spinors of NDE and the amplitude functions 
((nni)-(inHl) we have: 



= Fy = = tl^, (130) 



= -^^±^., Gy = -1±-1±, = (131) 



2i ' " 2i 

Substituting the solutions dHll) with spin j = 1/2 in ((T^ - ((nT| ) for 

arbitrary real a(r) and h{r) we obtain: 



i^x = -a(r)-; = -a(r)^; = -a(r)-, (132) 



= 0; Gj, = 0; = 6(r). (133) 
We rewrite again the conditions for the spatial phase functions: 

X F = 0; X (3 = 0. (134) 

These relations for solutions of kind H132|) - (|133() are satisfied only 
when: 

g = kor or g = kof{r), (135) 

and 

q = koz or q = kofiz), (136) 

where ko is the carrying wave number. The spatial phase functions of 
kind g = kQr and q = kQZ correspond to spectral limited pulses which 
satisfied additional relations: 
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AkAr = const. 



(137) 



The spatial phase functions of kind g = kof{r) and q = kQf{z) corre- 
spond to phase modulated pulses and for them the relations (|137() is 
not satisfied. The Poynting vector can be expressed by the amplitude 
functions of the electrical and magnet field: 

S = E{x, y, z, t) X H{x, y, z, t) = 
exp {i {W{t) - K{r))) F{x, y, z, t) x G{x, y, z, t) (138) 

where W and K are scalar phase functions. Substituting the solu- 
tions with spin j = 1/2 in above expression we find that: 

S = exp {i {W{t) - K{r))) ( -a{r)b{r)^; a(r)b(r)-; o] . (139) 

V r r J 

We see that the Poynting vector S is one circulation vector for 
solutions with spin j = 1/2 and its divergency is zero: 

V-5 = 0. (140) 

The relation (|14flj) determine that the energy flow through arbitrary 
closed surface around our vortex solutions with spin j = 1/2 is zero. 
The relation (|139|) show that flow of energy of our solutions circulate 
in X, y plane. Now we can generalize the above results for solutions 
with spin j = 1/2: The vortex solutions with spin j = 1/2 without 
external fields are immovable and electromagnetic energy oscillated in 
X, y plane. The electrical field oscillating spherically, in 'r' direction, 
while the magnet field oscillating in z direction. In the same way 
was calculated the Poynting vector for solutions with spin j = —1/2. 
For them we obtain that V • 5 7^ and we expect that they are not 
stable. One exact investigation of stability request investigation also 
the perturbation of the Poynting theorem and will discussed in one 
next paper. 

11 Conclusion 

In this paper we derive a set of Nonlinear Maxwell amplitude Equa- 
tions (NME) for nonlinear optical media with and without dispersion 
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of the electric and magnetic susceptibility. We have shown that in 
cases of linear and circularly polarized components of the electric and 
magnetic fields, the NME reduces to the Nonlinear Dirac system of 
equations (NDE). The equations are represented in a spinor form. 
Using the method of separation of variables, exact vortex solutions 
for both cases have been obtained. The optical vortex solutions ad- 
mit classical orbital momentum I = 1 and classical own momentum 
j = ±1/2. Here I would like to explain more abbout the differences 
between solutions, obtained by separation of variables of the usual 
linear Dirac equation with potential depending only on 'r' (for exam- 
ple hydrogen atom) and solutions of the NDE. In case of linear Dirac 
equations with potential to higher order spherical spinors correspond 
higher order of radial spherical Bessel functions. On the other hand 
for NDE to higher order of spherical spinors {i = 1,2, ..) correspond 
higher number of field components and higher value of localized en- 
ergy. For all radial solutions of NDE the zero spherical radial Bessel 
function ^lil^ are valid. The optical vortices in a media without dis- 
persion admit infinite energy integral. The energy integral of the vor- 
tex solutions is finite only in some special cases of paramagnetic media 
with suitable conditions on linear electric and magnet dispersion . Us- 
ing the Poynting vector for solutions with spin j = 1/2 we find that 
the energy flow through arbitrary closed surface around our vortex 
solutions is zero and the localized energy of our solutions circulate in 
X, y plane. Other important result is, that the vortex solutions with 
spin j = 1/2 without external fields are immovable. The electrical 
field in the vortices oscillating spherically, while the magnet field os- 
cillating in z direction. The initial investigations on stability of these 
solutions show the following: While the vortices with spin j = 1/2 are 
stable, the vortices with opposite spin (charge) j = —1/2 are not. 

All of the above results will be discussed later in relation with 
nonlinear field theory. However one direct approach to the plasma 
physics can be found. We consider the old idea presented by Gaponov 
and Miller in [^: Confinement of charged particles in potentials of 
high frequency electromagnetic fields. In our article are presented 
solutions inspired by this idea. The optical vortices are of such type 
of high frequency fields with zero intensity in the centrum. In a cold 
plasma it is possible to obtain trapping of ions in the field of vortices. 
Thus, conditions for ruled nuclear fusion can be reached. 
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